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' Generalized parton distributions (GPDs) are basic ingredients in the description of hard exclusive processes (see 
ref. and references therein) . Not only are they generalizations of the well known parton distributions from inclusive 
deep inelastic scattering (DIS), but being correlation functions they also incorporate non-trivial behavior of hadrons 
related to the nonperturbative regime of quantum chromodymanics (QCD). At present, apart from the preliminary 
studies of lattice QCD 0, one cannot calculate GPDs from first principles, so one has to resort to models or 
, parametrizations. 

fS| ' In fact, the perturbative approach to QCD is able to connect observables at different resolution scales, but the 
t realization of the complete project (i.e. to fully understand the consequences of the dynamics of quarks and gluons) 

■ requires the input of unknown nonperturbative matrix elements to provide absolute values for the observables at 
[ any scale. In the present paper we intend to apply a radiative parton model procedure which, starting from a low 

resolution scale Qo' been able to reproduce and predict (see, e.g., ref. Q and references therein) the main features 
' of the experimental deep inelastic structure functions at high momentum transfer. The procedure assumes that there 

Q-i, exists a scale where the short range (perturbative) part of the interaction is negligible and, therefore, the glue and sea 
^ are suppressed, and a long range (confining) part of the interaction produces a proton composed by (three) valence 

^ \ quarks, mainly Jaffe and Ross 6j proposed to ascribe the quark model calculations of matrix elements to that 

■ hadronic scale Qq. In this way, quark models summarizing a great deal of hadronic properties may substitute for 
low-energy parametrizations, while evolution to larger is dictated by perturbative QCD. 

In the following we study the nucleon's GPDs within specific hadron models and address the problem of evolving 
the input distributions to the experimental scale investigating the effects of different dynamical assumptions. In 
particular, we want to investigate both the quark core structure of the nucleon and its chiral properties. In fact, the 
new aspects of the GPDs with respect to the better known parton distributions are related to the socalled ERBL 
region where the presence of dynamical qq pairs, both in the nonperturbative and perturbative regimes, plays a crucial 
role. 

The importance of the chiral structure of nucleons is well established both experimentally and theoretically. The pion 
cloud associated with chiral symmetry breaking was first discussed in the DIS context by Feynman and Sullivan 1| . 
It leads to flavor symmetry violations in the sea-quark distribution of the nucleons naturally accounting for the 
excess of d (aiiti)q uarks over u (anti)quarks as observed experimentally through the violation of the Gottfried sum 
rule fTol[Tlllo.ll3| |. As discussed by Melnitchouk et al. the relatively large asymmetry found in these experiments 
implies the presence of nontrivial dynamics in the proton sea which does not have a perturbative QCD origin. In 
particular, a quantitative description of the entire region of the quark momentum fraction x covered by the experiments 
requires a delicate balance between several competing mechanisms. At larger x the dynamics of the pion cloud provides 
the bulk of the d — u asymmetry with DIS from the ttN component of the nucleon wave function, however also the 
ttA arises in the light-cone formulation of the meson-cloud model and it is of some importance in DIS too. 

Although the nucleon's nonperturbative antiquark sea cannot be attributed entirely to its virtual meson cloud |l5j . 
the role of mesons in DIS is of primary importance, and the idea was developed further giving origin to the meson- 
cloud model (for a review of early work, see refs. |16L [l7l and references therein). The connection between this 
model and the chiral properties of QCD was established by investigating the nonanalytic behavior of the d — u 
distribution (see, however, [l^). The meson-cloud picture is also suggested by QCD in the limit of large numbers 
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of colors iVc, where it becomes equivalent to an effective theory of mesons, in which baryons appear as solitons, i.e. 
classical solutions characterized by a mean meson field |2l| . A realization of this idea is achieved in the chiral quark- 
soliton model, where the effective action is derived from the instanton vacuum of QCD, thus providing a microscopic 
mechanism for the dynamical breaking of chiral symmetry p^ . Flavor asymmetry of the antiquark distributions 
arises in this model of the nucleon as a sub-leading effect in the limit of large Nc [23 • 

An alternative approach to investigate the role of gg-pairs in DIS and to access the ERBL region, is considering the 
constituent quarks as complex systems \24\ . Such a scheme has been recently developed in relation to a nonrelativistic 
constituent quark model, both for parton distributions J5] and GPDs 26]. 

In the present work we will study the possibility of integrating meson-cloud model effects into the evolution of 
GPDs. To this end we will assume that GPDs can be written in terms of double distributions |23|, involving a given 
profile function and the forward parton distribution derived in some model. At the same time the model proposed 
by Melnitchouk et al. 0| will be adapted to show how the knowledge of the meson-cloud effects can be incorporated 
within a relativistic quark-valence approach to GPDs. The double distribution (DD) model will be briefly recalled 
in sect.m Here the input parton distribution is discussed both in terms of the pure valence contribution derived in 
light-front relativistic quark models fsect. HT^ and in the meson-cloud model fsect. HTB)! . Matching sea, gluons and 
valence-parton distributions in QCD evolution of the obtained GPDs is then briefly described in sect. IIIII and the 
results are discussed in sect. IIVI Some conclusions are drawn in the final section. 



II. MODELING GPDS WITH DOUBLE DISTRIBUTIONS 



In the following we shall concentrate our attention on the chiral even (helicity conserving) distribution H'^ix, ^, Q^, t) 
for partons of flavor q at the hadronic scale where the models we are going to discuss are assumed to be valid to 
evaluate the twist-two amplitude. Such amplitude occurs, for example, in deeply virtual Compton scattering where 
a lepton exchanges a virtual photon of momentum q^^ with a nucleon of momentum P^, producing a real photon of 
momentum q'^ and a recoil nucleon of momentum P'^ . Then = —Q^Qfj, is the space- like virtuality that defines the 
scale of the process. In a symmetric frame of reference where q^ and the average nucleon momentum P^ = ^{P^+P'^) 

are collinear along the z axis and in opposite directions, the quark light-cone momentum is = xP^ , the invariant 
momentum square is t = = (P'^ — P^)^, and the skewedness ^ describes the longitudinal change of the nucleon 

momentum, 2^ = — . 

For sake of simplicity we follow the common notation which explicitly indicates three variables only {x, ^, t) instead 
of four (af, ^, We shall come back to the definition of when we discuss the values of the hadronic scale Qq 

and the QCD evolution of H'> in Q^. 

We also introduce non-singlet (valence) and singlet quark distributions, 

H^'^ix, t) EE ^ [m{x, t) + H'^i-x, e, t)] = t), (1) 
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tt) = Y. [^'(5^' ^' *) - ^'(-^' ^' *)] = -H^'i-x, t), (2) 

Q 

respectively. Besides being symmetric or antisymmetric in x, they are also symmetric under ^ ^ — ^ due to the 
polynomiality property j28| |. 

The analogous GPD for gluons is symmetric in x, i.e. 

H^{x,tt)^H\-x,^,t), (3) 

and reduces to the gluon density g{x) in the forward limit (x — > x) 

H3{x,0,0) = xg{x), x>0. (4) 

There are two distinct regions: the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) region, \x\ > ^, and the 
Efremov-Radyushkin-Brodsky-Lepage (ERBL) region, \x\ < The naming derives from the fact that the GPD 
evolution equations in the region \x\ > ^ {\x\ < £,) reduce to the familiar DGLAP (ERBL) equations in the limit ^ = 

(? = 1)-. 

The singlet and gluon distributions mix under evolution, while the non-singlet distribution does not. 
GPDs depend on the invariant momentum transfer t. In particular, the first moment inx oi H''{x, ^, t) is independent 
of ^ and related to the Dirac form factor of the proton. Thus in phenomenological constructions of GPDs it has been 
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found convenient to assume a factorized t dependence determined by some form factors. This simplifies the QCD 
evolution considerably because in this way the t dependences of quarks and giuons (which mix under evolution) are 
not modified during evolution. 

The t-independent part H'^{x,£^) = H'^(x,£,,t = 0) is parametrized by a two-component form p^ : 



where 



/•I ri-\0\ 

Hdd[x,0^ dp daS{x- l3-a^)F'}{P,a). (6) 

J-i J-i+\p\ 

The D-term contribution in eq. (0) completes the parametrization of GPDs, restoring the correct polynomiality 
properties of GPDs |2^. It has a support only for \x\ < ^, so that it is invisible in the forward limit. The 
D-term contributes to the singlet-quark and gluon distributions and not to the non-singlet distribution. Its effect 
under evolution is at the level of a few percent and in the following it will be disregarded in the input GPDs. 
According to Radyushkin's suggestion |23|, the DDs entering eq. © are written as 

F%(3,a) = h{f3,a)q{(3), (7) 

where q(0 ) is the forward quark distribution (for the flavor q), and the profile function h{P,a) is parametrized 
as [27im| 

h(3 a) = ^(^^+^) [{i-mr-^']' 

' ' 22b+ir2(fe + i) (1 - |;3|)2''+i • ^ ^ 

In eq. ©, the parameter h determines the width of the profile function /i(/3, a) and characterizes the strength of the 
f dependence of the GPDs. It is a free parameter for the valence ihvai) and sea {hsea) contributions to GPDs, which 
can be used in such an approach as fit parameters in the extraction of GPDs from hard electroproduction observables. 
The favoured choice is b^ai = ^sea = 1-0, corresponding to maximum skewedness. With a similar assumption adopted 
for the gluon distribution one defines hgiuon — 2. The limiting case b ^ co gives h{l3, a) 5{a) h(0) and corresponds 
to the ^-independent ansatz for the GPD, i.e. H'i{x,C) ^ H''{x,S, = 0) = q{x), as used in refs. '32,'33'|. 

In order to explicitly calculate H''{x,£,) in eq. © knowledge of the parton distribution q{x) is needed. In the two 
following subsections details are given about the derivation of q{x). 



A. Parton distributions and light-front relativistic quark models 

Following the lines of ref. 34] in two recent papers [HE HBl we discussed a method to evaluate GPDs within light- 
front constituent quark models (CQMs) at the scale dominated by valence (constituent) degrees of freedom. The 
comparison of these calculations with predictions in the chiral soliton model and the MIT bag model, as well as the 
consistency with lattice results for the first moments of GPDs showed that all the phenomelogy for large x and small 
t could be studied within the assumed relativistic CQM. As a drawback of such an approach, the calculation was 
restricted to the region x > ^ and the generation of qq contributions could have a perturbative origin only. In contrast, 
within the DD-based model both ERBL and DGLAP regions of quark GPDs are populated with any input quark 
distribution q{x) (with or without sea contribution). In this subsection we only consider valence quarks leaving for 
the next subsection the case includin g th e sea. 

According to the approach of refs. |35ll36| the parton distribution in relativistic light-front CQM takes this simple 
form: 

q{x,lxl)^Y.^r,r,j ildhsij^kA |*L^l({fc.;A„T.})P, (9) 



where fc^ = (A:° -I- fc^)/\/2 is the quark light-cone momentum, and A/q = X^i y^-i +"^1 the free mass for the 

three-quark system. ^^^{{ki]\i,Ti\) is the canonical wave function of the nuclcon in the instant form; under the 
assumption that only valence quarks are active, it is obtained by solving an eigenvalue equation for the mass operator 
within relativistic CQMs. 
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In the following we will discuss results based on the mass operator for the hypercentral CQM [37|, i.e. 

3 



M = ^ + mf - - + y, (10) 



with ki — 0, and rrii being the constituent quark masses, y = y fP + is the radius of the hypersphere in six 

dimensions and p and A are Jacobi coordinates. For a discussion of the model see refs. [3^l38l |. 

The distribution (|5J) automatically fulfills the support condition and satisfies the (particle) baryon number and 
momentum sum rules at the hadronic scale /Iq where the valence contribution dominates the twist-two response: 

dxq{x,fil) = Nq, (11) 

with Nq being the number of valence quarks of flavor q, 

dx x[u{x, fil) + d{x, fil)] ~ 1, (12) 

and u{x, /ip) = uv{x, /Iq) and d{x, /Iq) = dv{x, /Iq), the up and down valence-quark distributions. 

B. Parton distributions and the meson-cloud model 

Let us now introduce the meson-cloud model to incorporate qq contributions into the valence-quark model of the 
parton distribution discussed in the previous section. 

The basic hypothesis of the meson-cloud model is that the physical nucleon state can be expanded (in the infinite 
momentum frame (IMF) and in the one-meson approximation) in a series involving bare nucleons and two-particle, 
meson-baryon states. Its wave function is written as the sum of meson-baryon Fock states 

\p) = Vz\p)bare + Y. f dy d^k^^ <j)B M (v , kl) \B {y , k^); M {1 ^ y , ~k^)) . (13) 

BM'' 

Here ipBAiiyjk'^) is the probability amplitude for the proton to fluctuate into a virtual baryon-meson BAI system 
with the baryon and meson having longitudinal momentum fractions y and 1 — y and transverse momenta k± and 
— fc_L, respectively. Z is the wave function renormalization constant and is equal to the probability to find the bare 
proton in the physical proton. 

The lowest lying fluctuations for the proton which we include in our calculation are 

p{uud) — > n{udd) TT^ (ud) , 



p{uud) — !■ p{uud) 71*^ i —^[dd 



uu 



V2' 

p{uud) IS.^ [uud) t:'^ (^—j^dd ~ uu]^ , 
p{uud) A° (udd) 77+ (ud) , 

p{uud) A'^^{uuu)Tr'{ud). (14) 

In DIS the virtual photon can hit either the bare proton p or one of the constituents of the higher Fock states. In the 
IMF, where the constituents of the target can be regarded as free during the interaction time, the contribution of the 
higher Fock states to the quark distribution of the physical proton, can be written as the convolution 



^'ip{^) = Yl I — fMB/p{y)qM (-) + [ — fBM/p{y)qB 
MB L-'^y w / Jx y 



(15) 



where the sphtting functions fMB/p{y) ^nd fsM/piy) ^ire related to the probabiHty amphtudes (pBM by 

fBM/p{y) = /ms/p(i -y)= dkl mniy, kl) . (16) 

•^0 X,X' ' 
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Here ^1) probability amplitude for a hadron with given positive helicity to be in a Fock state consisting 

of a baryon with helicity A and a meson with helicity A' (391] . It can be calculated by using time-ordered perturbation 
theory in the IMF. The quark distributions in a physical proton are then given by 

qp{x)^Zqf^'{x)+5qp{x), (17) 
where q^^"^ are the bare quark distributions and the renormalization constant Z is given by 

Z=l-Y, dyfMB/p{y)- (18) 



MB 

AAA' 







The amplitudes 4>B]\iiy^ ^1) ma-Y be expressed in the following form 



'I^BMiy^kD = J- ^ , ,,, GHBM{y,kl)Vj'iiAy,kl), (i9) 



where mn is the ph ysic al mass of the fluctuating hadron (in present case the proton, but the approach can be 
generaHzed (e.g. ref. [33)), and 

Ml,, ^ '±±!< + kl+^<i (20) 

is the invariant mass of the meson-baryon fluctuation. ViMpiyi k±) describes the vertex and contains the spin- 

tendence of the amplitude. The exact form of the Vj)^jp{y, k'j_) can be found for various transitions in refs. 0, 
Because of the extended nature of the vertices one has to introduce phenomenological vertex form factors, 
GHBAiiy, ^1), which parametrize the unknown dynamics at the vertices. We use the popular parametrization 

In order to calculate the meson-cloud corrections to the quark distributions we have to specify the coupling con- 
stants entering l^A'^'p(y, fc^) and the cut-off parameters Asm- We use the numerical values as given by 0, l39l |. 
i.e. dNNTr/^'^ — I'^.G and S^att/^^ = 11-08 GeV~^. The coupUngs of a given type of fluctuation with different 
isospin components are related by isospin Clebsch-Gordon coefficients, gpmr+ = "V^SppTrOj gpAOTr+ = ~5pA+7ro/V2 = 
9pA++-!T- / VS, with gNNTT — QppTr" and qnAtt — QpA+^-K-- The violation of the Gottfried sum rule and flavor symmetry 
puts also constraints on the magnitude of the cut-off parameters. The values Amb = 1-0 GeV and A^/s = 1.3 GeV for 
the TriV and ttA components, respectively, give contributions to the u and d which are consistent with the requirement 
that the meson-cloud component of the sea quark contribution cannot be larger than the measured sea quark and 
also with flavor symmetry violation . With this choice of the parameters the probabilities of the fluctuations are 
given by Pjv^/p = 13%, Pat^/p = 11%. 

In the following we will assume that at the lowest hadronic scale the bare nucleon is described by the relativistic 
quark model wave function formulated within the light-front dynamics and, as a consequence, that only valence 
partons will contribute to the partonic content of the bare nucleon |35l l37j . The full (nonperturbative) antiquark 
content will be generated by the meson-cloud mechanism described by eq. H15|) . The partonic content of the A and the 
pion will be consistently evaluated within the same scheme assuming light-front dynamics and valence contributions 
only. Within these approximations the meson-cloud corrections at the hadronic scale /ig read 

q,{x)^Zq'';^-{x) + [/^^/^(l_y) + /^,/p(l-2;)]gV^') 

J X y \y / 

+ [fN./piy] + /AVp(y)] 9A , (22) 

where qp = (u ~ uy + u,d — dy + d) include both valence and sea contribution coming from the meson-baryon 
fluctuations, while 9^°'''^ = {u'^'^'^ , dy^^) include the valence contribution only as discussed in section III Al The 
conservation of both momentum and baryon number sum rules is guaranteed by the correct formulation of meson- 
cloud approach, in particular by the momentum conservation due to the symmetry fsM/piv) ^ fMB/pi^ ^ y) in eq. 
(|15|l and by the renormalization factor Z of eqs. (|13|) . H17|) and (|18|l . 
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The nucleon 

In order to model the partonic content at the scale /i§ for the nucleon, the A and the pion, we make use of the light- 
front approach discussed in sect. Ill K\ and calculate the diagonal component of the generalized parton distributions, 
i.e. the inclusive parton distributions by means of 

(23) 

where g(a;,/iQ) is given by eq. 10. 
The A 

The calculation of the cloud contribution involves the explicit form for the parton distributions i^) of the A (see 
eq. (|22() '): we use the results of the relativistic model for the nucleon and the isospin symmetries: 



(24) 



UAo{x,fil) = lup{x,fil), dAo{x,fil) = 2dp{x,fil), 
UA-{x,fil) = 0, dA-(x,/^g) = 3dp{x,fil). 

The pion 

The canonical wave function of the pion is taken from ref. and reads 

^2; Ml, M2) = ^3/!02 (1M15M2IOO) exp (-fcV(2/3')), (25) 

with k = ki = -fe, X = xi = k+/M^ X2 = I - x, M§ = {kj_ + ml)/x + {k^ + m2)/(l - x), and (3 = 0.3659 GeV. 
The choice of the model from ref. |4j| is consistent with the hypercentral CQM we adopt for the nucleon, in fact 
the central potential between the two constituent quarks is described as a linear confining term plus Coulomb-like 
interaction. The canonical expression H25|) represents a variational solution to the mass equation. 
The light-front parton distribution of the tt"*" is given by 

{x) = v^{x, fil) = ^ Sr^r^ Y{dh6 i^kA S ix - 
3=1 i=l \i=l / \ 

x|vI/W({fe,;Aa)f |(ir4-r,.|00)f . (26) 



Isospin symmetry imposes — dy^ — Uy — dy — u^(a;,/iQ), while, due to the model restrictions, the pion sea at 



the hadronic scale vanishes: u'^^ = d^^ = = d'^ = 0. 



III. MATCHING SEA, GLUONS AND VALENCE-PARTON DISTRIBUTIONS IN QCD EVOLUTION 

In order to extract the parton distributions of the nucleon including the sea (cloud) contributions we need to match 
sea, valence and gluons within the radiative parton model and to identify the matching scale Qq consistent with QCD 
evolution equation • 

We assume that the gluon distribution has a valence-like form and reads 

G{x, Ql) = ^ [uv{x, Ql) + dv{x, Ql)] , (27) 

where Ng represents the number of gluons. Since / dx [uy (x, Qq) + dv{x, Qq)] =2 + 1 because of baryon number 
conservation, we have 

j dxG{x,Ql)^^ j dx[uv{x,Ql) + dv{x,Ql)]=Ng. (28) 
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Following Gliick et al. |42|, we fix Ng = 2, the minimum number of gluons one needs to build a color singlet state. 
In spite of the simplicity of the assumption in eq. H27|l . the analysis of refs. [4^143^ shows that the crucial ingredient 
in the formulation of the radiative parton model is the amount of gluon momentum more than its actual form. Of 
course, to be really predictive one has to fit the form of the gluon distribution in a quite precise way, an accuracy 
which goes beyond the aim of the present study. 

The total momentum carried at the scale Qq (the scale of the physical proton, which will be larger than the scale 
of the bare proton fi^ made up of three valence only) must fulfill the requirement 



and consequently 



dxx [G{x,Ql)+uv{x,Ql)+dv{x,Ql)+2 [u{x , Ql) + d{x , Ql)]] = 1, 



dxx [uvix, Ql) + dv{x, Ql)] 



(29) 



1 + 
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dx2x[u{x,Ql) + dix,Ql)] 



0.52. 



(30) 



In conclusion, by using the notation (/)„ = J dxx^ ^ f{x) for the moments, we have {uy + dy)2 = 0.52, {G)2 = 
0.35, {2{u + d^)2 = 0.13, consistent with a next-to-leading order (NLO) evolution of the moments (in the DIS 
renormalization scheme) with Ql = 0.27 GeV^ and Anlo = 0.248 GeV [HII^. 

The actual values of Ql are fixed evolving back unpolarized data, until the valence distribution matches the 
required momentum {uy -\- dv)2 = 0.52. The procedure makes use of the valence contribution only, and it does not 
depend on the renormalization scheme. The value of Anlq is suggested by the analysis of Gliick et al. H^j, the 
coupling Q;s(Qo)|nlo obtained evolving back the valence distribution as previously discussed, and Ql is found by 
solving numerically the NLO transcendental equation 



In 



Ql 



An 

/3o cts 



PI 



In 



47r 



Po O-s 



Pi 
PI 



0, 



which assumes the more familiar expression 

«.(Q2) 



47r 



/3oln(QVALo) 



^ PI 



Po Inln(QVAgjLo) 
In(QVAko) 



(31) 



(32) 



only in the limit Q^ 3> A^^q (an interesting discussion about the effects of the approximation (|32|l can be found in 
ref. 0)- 

The hadronic scale, /Iq, consistent with the presence of valence degrees of freedom only (as discussed in sect. Ill A)l . 
is much lower and consistent with the constituent quark mass, its actual value being /Zg = 0.094 GeV^. The NLO 
evolution of the unpolarized distributions is performed following the solution of the renormalization group equation 
in terms of moments, i.e. {f{Q'^))n = Jq dx f{x,Q^) x"~^ , and involves kernels which have been computed up to 
next- to-leading order (NLO) in perturbation theory Since, in our case, the starting points for the evolution 

{^J^l,Ql) are rather low, the form of the equations must guarantee complete symmetry for the evolution from /ig to 
Q^ 3> fJ,l and back avoiding additional approximations associated with Taylor expansions and not with the genuine 
perturbative QCD expansion j43l l44j . 

In Fig. n we compare our results for the d — u difference with the data from the E866 experiment [l^. A QCD 
evolution with an SU(6) symmetric input introduces a very small amount of asymmetric sea at the experimental scale. 
Within our approach such an assumption corresponds to the perturbative contribution coming from three valence 
quark distributions at the lowest hadronic scale as predicted by the relativistic model of eq. ^ ,35j . The presence 
of the nonperturbative sea introduced by means of the meson-cloud model discussed in the present section improves 
quite substantially the comparison with the experimental data. In particular let us note: i) the important role played 
also by perturbative evolution in the region a: < 0.1 once the nonperturbative sea and gluon content is introduced at 
the hadronic scale Qg; ii) the satisfactory result in spite of our simple model. All evolutions have been performed 
according to ref. |43| 



IV. RESULTS AND DISCUSSION 



Results are presented in this section according to the model based on DDs as described in sect.m In particular, we 
will take: i) forward parton distributions with only valence quarks f section [ll A|l at input scale ^l = 0.094 GeV^ for 



8 



the DDs, and ii) parton distributions generated including meson-cloud corrections as in eq. (|22|l and matched with 
the gluon distribution at an input scale Qq = 0.27 GeV^, as described in section Hi Bl In both cases our input GPDs 
are continuous functions all over the whole range x G [—1, 1]. The t dependence is dropped from the very beginning 
and could be reintroduced in the final results by an appropriate ^-dependent factor. The D-tcrm is omitted as well. 

In addition to the results at the hadronic scale we will discuss the evolution of GPDs at higher scale. The QCD 
evolution was numerically performed to NLO accuracy (in the MS scheme), according to a modified version of the 
code of ref. ^| (see also ^Blp?!). The modifications we introduced in the main code are basically due to the need 
of a complete NLO evolution which makes use of the correct NLO transcendental equation H31|) . The original code 
makes use of the simplified expression (|32() largely unsatisfactory for our pourposes. 

In fig. 121 results are shown for the singlet quark, non-singlet quark and gluon GPDs obtained with no initial gluons 
and an input q{(3) in eq. ijTJl simply given by the bare (valence) parton distribution ^^'"''^(a:), eq. ©, derived from 
the hypercentral CQM. The model already gives a nonvanishing contribution to quark GPDs in the ERBL region at 
the hadronic scale /Xq without introducing discontinuities a.t x = ^ and with a weak ^ dependence. In particular, the 
absence of the sea contribution gives = H^^ at x > ^. After evolution up to = 5 GcV^ GPDs are almost 
confined into the ERBL region with a significant ^ dependence. 

This behavior is in agreement with previous studies of the GPD evolution showing that as the resolution scale 
increases the distributions are swept from the DGLAP domain to lie fully within the ERBL region ji^ . Functions 
with support entirely in the time-like ERBL region \x\ < ^ are never pushed out of the ERBL domain. In fact, the 
evolution in the ERBL region depends on the DGLAP region, whereas the DGLAP evolution is independent of the 
ERBL region. 

The same GPDs obtained when the input is implemented with the sea according to the meson-cloud model, i.e. 
qp{x) of eq. (|22|l . are shown in fig. 01 Now, the gluon GPD does no more vanish at the initial hadronic scale, that 
has to be redefined at Qq = 0.27 GeV^ according to the conclusion of sect. IIIII After evolution the qualitative result 
is similar to the case without the sea in spite of a more evident ^ dependence of the singlet quark GPD at the input 
hadronic scale. However, the shape of the non-singlet quark GPD around x = is sensitive to the input, and the ^ 
dependence of the GPDs size is in general less pronounced. 

Results have been presented for QCD evolution up to = 5 GeV^. This is already a value where GPDs have 
reached a stable configuration with respect to their asymptotic shape. In fact, the largest effects of evolution modify 
the input GPDs within the first few GeV^ in the evolution, as can be seen in figs. 14161 where the singlet quark, 
non-singlet quark and gluon distributions are plotted at £ = 0.2 as a function of Q^. A similar behaviour has been 
also found for the valence u GPD in the model of Ref. |2g . 

The results discussed till now have been obtained with the usual choice b^ai = bgea = 1 and bgiuon = 2 (compare 
eq. © and the discussion of sect.^J. It is known 113,1^ that in the DD-based model GPDs around x ~ ^ depend on 
forward densities at x <C ^, with a special sensitivity to the sea quark density. This is particularly critical when using 
q(x) and xg(x) taken from the global parton analyses [s^ because their singular behavior at a; = is responsible for a 
substantial and unrealistic enhancement of the quark singlet GPD relative to parton density in the DGLAP region at 
X ~ ^. In the model adopted here this singularity does not occur. Nevertheless, introducing the sea has nonnegligible 
effects on the size of the ERBL response, especially for the gluon distribution, as can be appreciated by looking at 
fig- [3 where the same curves for ^ — 0.1 from figs. [3 and Oar e directly compared. 

The value of the parameter 6 in the profile function h{j3, a), eq. JSJ, determines its width and has been shown to be 
inversely related to the enhancement of the singlet quark GPD 47, 50, 51] . In principle, one could choose h varying 
from the adopted values up to infinity, with a Dirac-delta-like width and quark GPDs reducing to parton densities. 
The h dependence of GPDs in the DD-based model including the sea contribution according to the meson-cloud model 
is illustrated in figs. |H1 and |51 for ^ = 0.1 and ^ = 0.5, respectively. Only the ERBL region is affected by the choice of b. 
At the input hadronic scale the h dependence is strong, increases with ^ and affects all GPDs, not only the non-singlet 
quark GPD. NLO evolution greatly reduces such a h dependence, and a similar behaviour is found for the evolution 
of GPDs also at LO. 

As a final remark let us discuss the validity of our evolution procedure. As already mentioned the evolution has 
been performed by using a code due to Freund and McDermott and modified to restore the correct NLO coupling 
constant thorough the transcendental equation (|31|l (see discussion in the previous section). The modification plays 
a crucial role when the evolution starts from low hadronic input as in the present investigation. With the caveat that 
perturbative stability can be tested only at the level of physical cross sections or observables like the unpolarized 
structure function F2 (see, for example, Ref. [s^), we show in Fig.^|the results of LO and NLO evolution for singlet, 
nonsinglet and gluon distributions up to = 5 GeV^ starting from a hadronic scale as low as Qq — 0.27 GeV^ 
(where sea contributions are included). In the the singlet and nonsinglet sectors, the results at LO and at NLO are 
quite close, showing that the evolution is clearly under control and converging. Similar results hold also in the case 
of the evolution of GPDs from the lower scale of /Xq = 0.094 GeV^, without the sea contribution. Only the gluon 
distribution shows larger discrepancy between LO and NLO evolution. This is not a drawback of our model as can be 
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seen from analogous results of Ref. 6u\ which make use of the standard (GRV98) |52'| parametrization for the diagonal 
partonic input at Qq = 0.36 GeV^, comparable with our input scale. In fact, it is well known that also in the case of 
parton distribution functions generated from low scale parametrizations, NLO gluons are affected by renormalization 
scale dependence, a problem not yet addressed for GPDs and deserving further investigation. 



V. CONCLUSIONS 



Different inputs at the hadronic scale have been considered in the QCD evolution of GPDs to study sensitivity of the 
results to the nonperturbative nature of the low-scale hadronic structure. In particular, the meson-cloud model was 
assumed to include sea quarks in the partonic content at the hadronic scale. Matching the sea and gluon distributions 
with the valence-quark distributions derived in a relativistic CQM, one starts evolution with continuous functions all 
over the range —1 <x< 1. From the results presented here we can draw the following conclusions. 

As already noticed in previous analyses |43, |43, |5J| , evolution in the DGLAP region is not much sensitive to the 
input. The reason is twofold. At large x the valence contribution dominates at the input hadronic scale and, as the 
scale increases, the distributions are swept from the DGLAP domain to lie entirely in the ERBL region. In addition, 
evolution never pushes the input distributions from the ERBL to the DGLAP region which then evolves independently 
of the ERBL input. In contrast, the ERBL region is rather sensitive to the input including or not the sea. 

Due to the mixing between singlet and gluon distributions in the evolution equations, even with a vanishing gluon 
distribution at the hadronic scale one obtains a significant gluon GPD in the ERBL region after evolution, peaked 
around x = 0. The peak height depends on the input and is higher when the sea is included. 

The present results focus on the ERBL region as the most interesting one to look at the nonperturbative effects 
surviving evolution. This is suggesting that one has to investigate suitable processes under appropriate kinematic 
conditions to study such effects. An analysis of possible observables is in progress. 
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o E866 exp 




FIG. 1: d — u distribution: at the scale Qq = 0.27 (GeV/c)^ (dotted curve), at the (average) scale = 54 (GeV/c)^ of the 
E866 experiment (solid curve), at the comparison scale of = 5 (GeV/c)'^ (dot-dashed). The result (multiplied by a factor 
10 for convenience) of an evolution starting from the lowest hadronic scale where the nonperturbative sea content is neglected 
is also shown (dashed). Data from the E866 experiment [l^ . 
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FIG. 2: Singlet quark (upper panels), non-singlet quark (middle panels), and gluon (lower panels) GPDs using input parton 
distributions of the hypercentral CQM. Full lines for ^ = 0.1, dashed lines for ^ = 0.2. Left panels refer to the input hadronic 
scale //o, right panels to the corresponding results obtained from NLO evolution up to = 5 GeV^. 
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FIG. 4: Singlet quark GPD at ^ = 0.2 as a function of x and obtained from NLO evolution of the GPDs using input parton 
distributions of the hypercentral CQM with the sea contribution at the initial scale Qq. 



FIG. 5: Same as in fig-BI but for the non-singlet quark GPD. 




FIG. 6: Same as in fig.ll but for the gluon GPD. 
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FIG. 7: Singlet quark (upper panels), non-singlct quark (middle panels), and gluon (lower panels) GPDs at ^ = 0.1 obtained 
from NLO evolution up to = 5 GeV'^. Full lines refer to input parton distributions of the hypercentral CQM at the lowest 
hadronic scale /liq = 0.094 GeV^, dashed lines refer to results including the sea contribution at the hadronic scale Qo = 0.27 
GeV^ 
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FIG. 8: Variation with the profile parameter b of the singlet quark (upper panels), non-singlet quark (middle panels), and 
gluon (lower panels) GPDs derived for ^ = 0.1 using input parton distributions of the hypercentral CQM and including the sea 
contribution. Left panels refer to the initial scale Qq, right panels to the corresponding results obtained from NLO evolution 
up to = 5 GeV^. Full lines for 6=1, dashed lines for b = 10, and dotted lines for & — > oo. 
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FIG. 10: Singlet quark (upper panels), non-singlet quark (middle panels), and gluon (lower panels) GPDs using input parton 
distributions of the hypercentral CQM with the sea contribution in the meson-cloud model, evolved at = 5 GeV^ and at 
^ = 0.2. The full lines are the results from NLO evolution, and the dashed lines corresponds to LO evolution. 



